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Abstract 


Silicon forms the backbone of modern electronic industry. The crystalline state of this 
element in its pure or impure form, is the most commonly used material. However other 
forms of silicon are also increasingly finding their use, and have already established a niche 
for themselves. Thus amorphous silicon passivat('d with hydrogen (a-Si;H) has emerged as a 
popular optoelectronic material with widespread application in solar photo- voltaic industry. 
The material has also revolutionized the growth of microelectronics with applications in 
LCD/TFT's and light sensors. But inspite of the increasing technological implications of this 
material, not much is known about its physics. The issues of thcuunal instability, sensitivity 
to visible radiation, affecting the properties of a-Si:H are yet to be unch'rstood. 

In this study we have focussed our attention on hydrogenated clusters of silicon. Clusters 
constitute an assortment (from 2 to 100) of silicon atoms in a.ssociation with hydrogen. They 
exhibit markedly different material propertic's from !)ulk silicon with potential application in 
cluster assembl(;d material. Mon.'ovcu' a study of these clusters is exp(H;tcd to throw light on 
the Staebler-Wronskii efffict exhibited by a-Si:H. Implications also exist in the realm of porous 
silicon another futuristic material with potential application in LEDS and semiconductor 
lasers. Amorphous silicon and porous silicon together has raised the hope of producing 
integrated opto-electrical circuits. 

Calculating the ground state geometry of clusters was the principal objective of this 
study. This constitutes in optimizing structures of clusters generated randomly, using genetic 
algorithm. The cohesive energy of cluster comprised of our objective function. A confined 
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cubic space with specified dimension was chosen (Cubic box with dimension 4 to tA) wherein 
the cluster evolves, giving the final ground state configuration. Thus for a n atom cluster 
the total number of variables to be optimized is 3n, consisting of x, y, z, coordinates for each 
atom. The lower and upper bounds of each variable is given as 0< Xj < L where L is the 
length of cube side. 

The cohesive energy of a cluster is calculated using the principles of Kon-orthogonal 
tight binding theor}-. According to this theory, cohesive energj' of a system comprises of two 
principal parts: a classical contribution, consisting of coulombic interaction between atoms 
and a quantum part involving electronic interaction in the atoms itself. This particular 
methodology was chosen since it is both accurate and also computationally efficient. It 
is especially true for the Si-H semiconductor system. In this work we have prcsenti'd the 
optimized structure's of a number of a-Si;Il clusters for (n = 1-2, H = 1-5) and Si^H 

(k = 3-8) varieties. 
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Chapter 1 


Introduction 


Hydrogenated amorphous silicon (a-Si:H) is increasingly becoming an important mate- 
rial, emerging a long way from its initial beginning as a potential candidate for solar cell 
applications [1], Since then interest in the material has grown up rapidly with its application 
extending into the field of large area microelectronics, sensors, detectors etc. The various 
advantages of amorphous Si-H and its alloys are discussc'd below 

• Short optical absorption length for visible light. 

• Good intrinsic and doped semiconductor properties. 

• Deposition easily at low cost over extremely large areas. The temperature of deposition is 
low and is of the order of 250 "C. 

• Monolithic integration of individual cells is simi)le which avoids costly cut and paste meth- 
ods. 

• Both silicon and hydrogen arc abundant material constituting more than 50 -60% of the 
earths crust. 

Due to the above mentioned qunlities a-Sidl has become the most used thin film material 
for photo-voltaic cell constructions and as contact image-sensors in scanning machines. A 
schematic diagram of a contact image-sensor is shown in Fig. 1.1 [2]. 

The ability of a-Si:H sensors to copy over large paper width, high resolution of imaging 
and simplicity in optical system design, has made them useful for document scanning. 
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a-Si;H Scan 
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Figure 1.1: Schematic of electrophotographic sensor (Source: [2]) 

Due to their high sensitivity to long wavelength visible light, high dark resistivity and 
strong resistance ro mechanical wear, a-Si;H alhjys have also been used as photo-receptors for 
electro-photography. Their first commercial use in ('l(.'ctroi)hotographic drums was done by 
Canon Inc. in 19S4. The smaller band gap of a-Si:H alloys compared to amorphous selenium 
allows a shorter discharge time, which implies higher copying spei'd. In electrophotographic 
processes the thickness of a-Si:H films required is of the order of 20//m. This reciuirement 
of thickness is an important issue to be considered since the deposition rates of high-quality 
a-Si:H films are below 0.4 nm/s. 

Amorphous Si:H has also found its utility in photodiode applications like ultraviolet 
light, X ray and charged particle detectors, as image sensors for neural network applications 
[3] etc. 

The advantages of producing uniform film coatings, high yield and consiflcra!)le field 
effect mobility offered by a-Si:H, has established it as an ideal material for large area elec- 
tronics applications. A widely used application of a-Si:H is in the fabrication of thin film 
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Figure 1.2: Cro.s.s-Section of an a-Si:H / color Display (Source: [4]) 


transistors (TFT) used in controlling liquid crystal displays (LCD) [4]. A schematic of such 
a device is given in Fig. 1.2. Presently this technology is heavily used for the manufacture 
of portable television and in inilitary/avionic api)lications. 

The fast progressing technokjgy of a-Si:hI is h<j\vever plaguc^d by its own characteris- 
tic problem. All amorphous Si:H alloys exhibit reversible df'gradations in many material 
characteristics on being subjected to light exposure, charge injection or thermal (pienching. 
Understanding the physics behind the.se phenomenon is therefore becoming increasingly im- 
portant to give rise to a more rcdiable technoh^g}'. In the next section we have outlined, the 
common techniques used for producing device quality a-Si:H. 


1.1 Production of a-Si:H Based Materials 


1.1 Production of a-Si:H Based Materials 

The different deposition techniques for a-Si;H are [5]. 

• Chemical vapor deposition technique which includes , PECVD, PACVD etc. 

• Reactive sputtering. 

However the chemical vapor deposition technique employing glow discharge is the most 
popular method. In this method the starting gas material used is generally silane (S'fRT^). 
This gas (also known as carrier gas) is introduced at a controlled rate in between two elec- 
trodes in a evacuated gas chamber, where it is subjected to glow discharge produced by 
initiating a direct current (DC) or alternating current (radio frequency or microwave). The 
discharge causes the gas to dis.sociate and ultimately get d(?posited on a substrate producing 
a-Si:H film. A simplified diagram of tin' glow di.scharge process is given in Fig. 1.3. 

1.2 Staebler-Wronskii Effect 

Staebler and Wronskii [C] re-ported that a-Si:H films on being (^xi)o.s(’d to visible radiation 
shows an alteration in their material proi)ertie.s. Thus the prope-rties of dark conductivity 
and photo-conductivity in a-Si:H show a diminishing trend on being expos(;d to light. This 
effect is stable at room temperature, but is reversed to its original state on annealing treat- 
ment at 180 °C for a time period of one to two hours. Treatment of the film by infrared 
light however fails to restore its initial state. Ev(tr since the discovery of this phenomenon 
lot of research work has been done to get a clear umhtrstanding of this effect. As of now no 
universally accepted model for tin- iindc-rlying nu-chauism has been found. However some of 
the widely accepted experimental facts in this regard are given below [7]. 

• This light-induced degradation is observed in all a-Si;H films regardless of its method of 
preparation and film quality. 

• The phenomenon is a bulk effect. I 

• Light soaking with photon energy larger than 1.1 ev, can all induce this effect. Thus 


1,2 Staehler-Wronskii Effect 


i 



Figure 1.3; Sduunatic of glow-discharge deposition process 

bombardment by high energy particles, (luenching from high temperatures, charge injection 
etc., have a similar effect on a-Si:H. 

• The effect decreases with tlu; decrea.se of sample temperature. 

• The defect can be removed by ihernial annc'aling at 100 to 200 °C for about one hour. 
Infrared radiarion has no effect on the jjroce.ss of annealing; visible' light however increases 
the annealing rate. 

• Neutral silicon dangling bonds formed during the creation of this metastable state are 
situated 0.8 ev above the valeiiGC band mobility edge. 

• The metastability is proposed to be an intrinsic property of a-Si:H alloys, since impurities 
below 1*10^® cm~^ have no detectable effects [8]. It is worthwhile mentioning here that 
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concentration of most impurities, like P, B, C, and N, can be reduced to levels lower than 
the concentration of light-induced defects without improving the stability of the material. 

• The phenomenon is not directly related to any intrinsic mechanical stresses generated dur- 
ing the film preparation stage [9]. 

Several models for the microscopic mechanism of metastability in a-Si:H alloys have been 
proposed. Some of the well known models are [7] 

1. Weak-bond-breaking model [10, 11, 12]. 

2. Charge-defect model [13, 14, 15]. 

3. Reliybridizcd two-site model [16]. 

Most of the models involve the presence of hydrogen in the material. ■ Experimental re- 
lationship between the diffusion coefficient of hydrogen and annealing rate has also been 
established. All this warrants a fundamental understanding about the molecular structures 
and internal bonding of a-Si:H, in particular. Furthermore, identifying the controlling reac- 
tions in the production of hydrogenated amorphous silicon involves a thorough knowledge 
about the various Si-Il clusters and their relative stal)ility. In this work we focus on Sijffim 
(n = 1,2 and m = 1,5) and SikH (k= 3,S) clusters. 

In an earlier study, based upon uon-orthogonal tight binding molecular dynamics model 
[17], Gupte and Pra.sad [IS] reported the ground state' structures and the vibrational spectra 
of these clusters using the simuhited annealing techniciiu'. Although a large number of 
structures wore obtained during that work, the reciuirernent of large CPU time often renders 
the usage of simulated annealing in such problems computationally prohibitive. Our main 
aim in this work is to replace simulated annealing by genetic algorithm (GA) [19] in order to 
calculate the ground state structures with better efficiency. The computational gains using 
GA can be extended further with the parallelization of our code. It is worth mentioning 
here that, GA is highh- suited for parallel processing unlike any other traditional or non- 
traditional optimization methods. In the next chapter the fundamentals of GA are discussed. 




Chapter 2 


Fundamentals of Genetic Algorithm 


Traditional optimization methods are derivative based. Starting from an initial guess, 
the search algorithm selects a new point based upon some derivative criteria, leading to the 
nearby optimum. Extensive details on traditional optimization methods are available else- 
where [20]. However in genetic algorithm the approach is entirely different in that a set of 
initial points, individual in GA parlance, is the starting point. These group of individuals 
population in GA terminology are generally subjected to three powerful genetic operators: 
reproduction, crossover and mutation [19] described later. The basic principle of genetic 
algorithm is based upon the quasi-Darwinian principle of survival of the fittest where an 
individuals claim to pass its clone to the next generation is evaluated in terms of its fitness, 
which is either the function itself or .some appropriate transformation of it. For example, 
in order to optimize a function 0(x) between the lower and upper bounds and x^\ the 
fitness function f(x) can be taken identical to ^(x) in a maximization case. On the other 
hand, in a minimization problem, the fitness function can be defined as 


/(x) = 


1 

1 -+■ <j){x) 


(2.1) 


Several other suitable transformation of the objective function can be taken depending upon 
the nature of problem to be solved. In GA the variables are generally coded as binary strings 
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using ’I’s and ’O’s. Thus if four bits are used to code the variable x, then the strings (0000) 
and (nil) would represent and respectively, since those are the strings having the 
minimum and maximum decoded values. To convert the binary variables to real, different 
transformation schemes are used. A commonly used rule is linear mapping which is stated 
as: 


X = 


x^-x^ 

2^-1 


(5) 


( 2 . 2 ) 


In the above equation S is the decoded value of the binary string s of length I, representing 
the variable x. If the j’th bit in the string is denoted by sj, then 


i-i 

5 = ^ 2^Sj where Sjc(0, 1) (2.3) 

j=0 

Thus the four-bit string (1111) Inus a decoded value of 16. All the variables representing a 
point arc mapped onto their binary equivalents which are concatenated together to form a 
single string whose fitness is then calculated. Using the initial population one then proceeds 
to create a mating pool consisting of individuals selected on the basis of their fitness values. 
The reproduction operator is used for this purpose. A detailed discussion on the various 
mechanisms of reproduction arc given below. 


2.1 Reproduction Operators 

Initially reproduction (or selection) operator is applied on a given population. The con- 
sequence of this operation is to select the higher fitness strings, which together create the 
mating pool. These strings in the mating pool arc then subjected to crossover and mutation 
processes before being passed into the next time step, generation in GA terminology. As in 
this real world, an individual in GA is subjected to the same phenomena of birth and death 
which can be expressed mathematically as [21] 
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h,t+i = Ik,t + h,t,b — h 


t,d 


( 2 . 4 ) 


Here I denotes the number of individuals with the same fitness fk, ’k’ being its fitness class 
identifier, while ’t’ is the time at which the class is considered. The subscript ’b’ denotes the 
individuals being born in that class while ’d’ symbolizes those dying. In a normal scheme 
(also known as overlapping population model) all the individuals existing in the present 
time step is considered to die, birth rate being the sole criteria determining the number of 
individuals existing in that class at the next generation. One of the more commonly used 
selection operator is the fitness proportionate selection method. Here the probability of 
selecting an individual of the the i’th class into the mating pool can be approximately given 
as 


Pi = 




( 2 . 5 ) 


Goldberg and Deb [21] have shown that for a non-overlapping population of constant size 
n the expected number of coj)ics of i’th class (among a total of k classes) in the following 
generation is given by 


PiMi = pA (2-C) 

Jt 

where ft — is the average function value at the current generation, n being the 

population size. Solving equation (2.0) we get [21] 


hPi,0 


( 2 . 7 ) 


In order to have an understanding about the convergence situation, Goldberg and Deb [21] 
undertook an analysis for the takeover time required by this selection method; takeover time 
being definc<l as the time re<iuired by an individual of the best group to acciuirc a ^ pro- 
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portion of the total population. Considering a polynomial objective function the calculated 
takeover time is found to be [21] 


= -(rrlnn — 1) (2.8) 

c being the exponent of the functions considered and is a constant. For an exponential 
function the takeover time is given as 

C = -ln?r (2.9) 

c 

Thus for exponential objective functions the proportionate distribution method converges 
in 0(ln n) time, a factor of ii faster than that required in the polynomial case. However a 
major problem in using this type of reproduction operator is premature convergence. During 
the initial time fitness variance in the population is high, with few high fitness strings while 
the remaining being of much lower fitness. On using this mechanism of reproduction the 
higher fitness strings propagate quickly resulting in a lack of diversity in the population [22]. 
Avoidance of such problem involves the usage of low selection pressure at the initial period 
when population variance is high and inc.rcasing the pressure at a later time period when 
the variance decreases. 

Sigma scaling technicpui incorporates this mechanism. Here an individual’s expected 
value is expressed as a function of its fitness, poimlation mean, and the poj)ulation standard 
deviation. An examples of sigma scaling is given l>elow [22] 


Ki:p Val(i,t) 


1.0 + if(T(t)y^() 

1.0 if (T(t) = 0 


( 2 . 10 ) 


where ExpVal(i,t) is the expected value of individual i at time t, f(i) being the fitness of i, 
/(t) the mean fitness of the population at time t, and (7(t) the standard deviation of the 
population fitness’ at time t. 




2.1 Reproduction Operators 


11 


Rank selection operator is another selection methodology where an individual is selected 
on the basis of their ranking. The takeover time for linear ranking scheme [21] is given as 

t* = -log{n — 1) (2.11) 

c 

In ranking selection the procedure involves sorting the entire population according to their 
individual fitnesses which is done on the basis of some assignment function. This kind of 
selection technique has been identified to be potentially time consuming, specially when the 
population size involved is large [22]. 

To avoid this disadvantage the tournament selection operator is the most appropriate 
method. It includes the advantage of rank selection operator in maintaining a uniform se- 
lection pressure and is computationally much more efficient being well suited for parallel 
processing. The implementation in\’olves selecting two individuals randomly from the given 
population which arc then allowed to compete against each other. A uniform random num- 
ber is g(inerated between 0 and 1 which is then compared to a parameter k; if the number 
is greater than k, the fitter individual is selected or otherwise. The tournament selection 
method is sometimes associated with Elitism where the fitter individual is selected preferen- 
tially. It maybe mentioned here that tournament selection with more than two competitors 
has also been devised. 

Goldberg and Deb [21] has .solved the takeover time for tournament selection which is 
expressed as 


r = -^[In n -+- ln(ln n)] (2.12) 

Ins 

where s is the tournament size. 

Of the different reproduction operators stochastic universal proportionate, stochastic 
remainder proportionate and tournament selection algorithms are the only ones with a 0(n) 
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time complexity and are the recommended selection operators [21]. In our work we have 
used the binary tournament selection opei'ator. The string with higher fitness was selected, 
the process being repeated ’n’ times ’n’ being the population size. In genetic algorithm, 
unlike in natarul biological systems, the population size is kept constant throughout. The 
essential idea of reproduction is to assign more copies of the string with higher fitness to the 
mating pool, consisting of selected individuals where the crossover and mutation operators 
are applied with a pre-assigned probability. Detailed discussion on the crossover mechanism 
is given below. 

2.2 Crossover 

Single-point crossover is the technique commonly used by GA researchers. Here two 
strings are randomly chosen from the mating pool. For example in a six-bit population: 

110100 

and 

001011 

could be two typically selected individuals. Next, a uniform random number between 0 and 
1 is generated. If the random number is less than a pro-assigned probability, the process of 
crossover is acti\-ated by selecting a random crossover site (j) and swapping the bits right to 
it. For example, if the random crossover site is after the third bit from the right hand side, 
such that 

no I 111 
001 1 on 

then the crossover process results in the offsprings 

noon 

and 


001111 
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In GA literature the above mentioned crossover procedure is known as single point crossover. 
Another variant of crossover is the two-point crossover technique. Here two crossing sites 
are chosen at random and the segments between them exchanged. Thus consider the two 
strings with their crossing site as given below 

11 I 01 I 11 
00 I 10 i 11 

After crossover the two new strings generated are 

111011 

000111 

As an extension of the two-point crossover technique multi-point techniques have also been 
developed. Uniform crossover method is a multi-point technique stretched to its extreme. 
Hero each bit of the individual is exchanged with a probability ’p’. The value of ’p’ is 
generally kept in the range 0.3 - 0.5. 

The selection of proper crossover techniciue is very important. However no simple so- 
lution exists to predict the proper choice of crossover technique, the issue being variously 
determined by the nature of function used, variable coding and other GA parameters. In a 
compar;itive analysis of different cro.ssover techniciues [22] the two point crossover method has 
been prcferrc'd. The single point crossover tecimiciue is said to be more prone to positional 
bias and endpoint effect. This occurs as a consequence to the fact that single-point crossover 
favors the formation of short schema ^ which is invariably accompanied by the preservation 
of unfavourable bits associated with short range schemata’s. Besides, in single point mech- 
anism the string ends are always exchangc'd during crossover. The experiments with Royal 
Road functions have also shown up these effects [22]. However it would be inappropriate to 
make a deterministic view about the choice of crossover techniques from empirical evidences 
obtained from test runs on a specific function. The uniform crossover method on the other 

^Schema (Schemata in plural) represents the commonality in bit positions among a group of strings. Thus 
in a three bit problem, the (10*) schema, displays the similarity among the (100) and (101) strings. For 
further reference check [19]. 
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hand is free from positional bias but at the same time is highly disruptive to schema forma- 
tion mechanism [22]. Theoretical studies on the search power of uniform crossover with its 
extension to more complicated cases are available in [23]. 

In the present study both single and two-point crossover methods have been applied. 
Uniform crossover technique was also applied, but was found to be inferior with respect to 
the performance of single and two-point methods. 

2.3 Mutation 

After crossover the mutation operator is applied. The idea is to fecilitate a local search 
by sparingly altering some of the bits. Mutation operator therefore examines each binary 
bit of an individual, and attempts to change it to its complement with a small pre-assigned 
probability. While examining any particular bit, a uniform random number is generated. 
If the random number is less than the mutation probability, the bit is changed to its com- 
plement. New strings created at the end of these processes arc evaluated for their fitness’. 
The mutation process described above is known as jump mutation] another variant called 
creep mutation is also used. In creep mutation each variable is examined in the real space, 
and is either increased or decreased by a small value with a probability of 0.5 which is then 
reconverted into the binary space. 

After the process of mutation is over the new strings obtained are evaluated for their 
fitness’ which completes a cycle of GA iteration. The same procedure is continued until a 
terminating criterion is satisfied. In the next few sections some advanced operators used in 
GA are discussed. 

2.4 Niching 

In biological terms niche can be explained as an organism’s job or role in an environment 
[19]. The importance of niching in the light of function optimization is discussed here. Let 
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us consider a function as shown in Fig. 2.1 which we want to optimize using a simple genetic 
algorithm. W'e begin with a random population on which the genetic operators of reproduc- 
tion, crossover and mutation is applied. It is likely that in such a case the population will 
tend to converge onto one of the peaks shown. This is principally due to the so called genetic 
drift resulting from stochastic errors caused by small population sizing [19]. Alternately a 
optimization problem can lead to a function exhibiting multi-modality with different peak 
heights. In many practical cases we are interested in identifying all the relevant peaks which 
is otherwise impossible to do using traditional optimizers, at least in a single run. The tech- 
nique of niching in genetic algorithm plays a major role in this regard. 

A number of techniques exists to implement the niching scheme in genetic algorithm. Gold- 



Figure 2.1: Function showing equal peaks (Source; [19]) 


berg and Richardson [24] utilized this concept by introducing the sharing function. Here 
the population is divided into different sub-populations according to the similarity of the 
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individuals [25]. This can be done either in the real space or even in the binary space itself. 
When the implementation is done on the real space it is called Phenotypic Sharing, and in 
binary space it is known as Genotypic Sharing. A parameter Cshare [24] is constructed, which 
determines the sharing existing between the strings. The power-law sharing function Sh(d) 
is then defined as [25] 


Sh(d) = 


^(^share'^°‘ if d < 

0*0 ^ ^ ^share 


(2.13) 


where d is defined as the distance between two individuals. Once the extent of sharing 
is determined an individual’s fitness is degraded accordingly. The working of the sharing 
principle depends chiefly on the parameter CTshare the choice of which is vital [25]. The 
parameter ashare depends on the number of peaks and the lower and upper bounds of the 
solution space [25]. 


2.4.1 Phenotypic Sharing 

In this type of sharing the distance (dij) between the strings can be calculated using 
any appropriate distance norm in the n-dimensional space [25]. The Eucledian distance in 
the n-dimensional space can also be adopted [25] and is given as 




Ifti 

M E (^k,i - Xkjf 

\k=l 


(2.14) 


where Xi,£, arc the parameter values in the real space for the fi’th individual. 

For estimating the parameter cr,.harc C'^ch niche is considered to be enclosed in a n- 
dimensional hyper-sphere of radius ashare such that each sphere encloses 4 of the volume of 
the solution space, q being the total number of peaks existing [25]. The radius of a hyper- 
sphere enclosing the entire space is then given as 
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r - - 


2tt 


,max ^kjTnin) 


(2.15) 


the volume of sphere being calculated as V = cr" where c is a constant [25]. Thus cXshare is 
then calculated as given below [25] 


^^share 


■cr ‘ 


O-shc 


vTETfe ytriax - ,min) 

__ 


(2.16) 


(2.17) 


2.4.2 Genotypic Sharing 


In genotypic sharing the genetic proximity between two individuals is defined by the 
number of different alleles in their chromosome [25]. Thus consider two strings Si and Sj 
both of string length /. Comparing the two given strings, if only one bit difference is allowed 
btdween them tlnm 'Ci is the total number of allowable strings. Thus for k allowable bit 
differences the total number of possilfie strings are As the total number of allowable 
strings in the binary space is 2^ so proportion of strings which arc k bits different is ^Ca:/ 2' of 
the given solution space. Thus the total proportion of strings with k bit differences or less is 
given by the exprc.ssion . Now assuming that the function considered has q peaks 

distributed uniformly in the entire variable domain, each niche should then correspond to an 
average of i of the total solution space. Thus k being the maximum allow'able bit difference, 
therefore 






t=0 


1 

Q 


(2.18) 
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2.4.3 Triangular Sharing Function 

In the present study Goldberg's triangular sharing function has been used (see Fig. 2.2). 
In this method the difference in distance between each parameter to be niched is found out 



CJjj =:||Xj - Xj li 

Figure 2.2; Triangular sharing function (Source; [19]) 

.for all the strings. This distance is generally taken in a normalized form. Thus for any string 
its share is obtained by summing up all the differences for the niched parameters. Eventually 
for a given string the total extent of sharing is determined by summing up the contribution 
of each individuals in the entire population. Thus any individual in proximity with the given 
string has a sharing function value close to one; the further away the string is, the sharing 
function value decreases correspondingly. On getting the total sharing function value the 
individuals fitness is degraded using the following scheme 
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2.5 MicroGA 

The simple genetic algorithm (SGA) discussed above uses only three basic GA operators, 
namely reproduction, crossover and mutation. This kind of GA is suitable for stationary 
function optimizations, but fails to optimize non-stationary functions where the function 
itself is evolving at a rate faster than SGA’s [26]. This kind of real-world problems include 
problems in aerospace, such as pursuit and evasion, on-line aircraft trajectory optimization, 
optimal control of aircraft in wind shear [26] etc. Krishnakumar [26] proposed a new scheme 
of genetic algorithm for application in such purposes and also in stationary optimization. 
The stepwise procedure for microGA implementation as proposed by [26] is given below 

1. A micro-population usually of size 5 is selected randomly including one good string from 
the previously conducted search. 

2. The individual fitnesses are evaluated and the best string is identified. The best individ- 
ual is then passed into the next generation (elitist strategy). This procedure guarantees the 
preservation of good schema. 

3. The remaining four strings are chosen using a tournament selection strategy. Here a de- 
terministic tournament selection method [19] is used since the law of averages does not hold 
good for small population sizes as u.sed in microGA. Simultaneously precautions .should be 
taken to prevent mating of same strings for the next generation. 

4. Crossover is then applied on the selected individuals. A cro.ssover rate of pc = 1 is to be 
used. This imi)lies that each and (!very string is subjected to crossover action. Using a high 
crossover rate is expected to promote schema processing; the mutation rate is kept at zero 
as enough diversity is introduced into the population through the cros.sover action. 

5. The population is then checked for convergence ( either genotypic or phenotypic conver- 
gence). If convergence is achieved we then proceed to step 1. 

6. Else to step 2. 
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2.6 Inversion 

In nature the functionality of a chromosome is independent of its gene positions. This 
is however not the case in genetic algorithm, where bit positions comes to play an important 
role. For enabling GA’s to mimic nature the inversion operator was developed. Here each 
gene in a chromosome is attached to an identifier giving the locus of that gene. Hence the 
identity of a chromosome can be preserved independent of the position of its genes. Thus 
in ordinary binary GA a four bit string 1101 is stored under inversion conditions as (1,1), 
(1,2), (0,3), (1,4). Under such a situation the position of a gone is no longer necessary to be 
fixed. The inversion operator in genetic algorithm is implemented by fixing two bit positions 
on the chromosome randomly and then reversing the scciuence of genes in between these two 
positions. Thus consider that the two bit positions selected for the 4-bit string shown above 
are situated in between the first and second and second and third bit positions respectively. 
After inversion the new string is given as (1,1), (0,3), (1,2), (1,4). The genetic operators 
of crossover and mutation are then allowed to operate on this string. It may be noted 
that the normal crossover operator under inversion conditions is redundant since this would 
entail the possibility of same genes being repeated in the new chromosome. Thus starting 
with two strings (1.1), (0,3), (1,2), (1,4) and (1,0), (1,2), (0.3), (1,4) with the cros.sover 
point chosen in between the 2'“^ and 3’’'^ bit positions crossover yields’ two offsprings (1,1), 
(0,3), (0,3), (1,4) and (1,0). (1,2), (1,2), (1,4). Thus the same genes gets repeated in the 
chromosome making the cros.sover redtmdant. To avoid this situation certain strategies have 
been proposed [19, 22]. One of the strategy is to allow crossover between strings which have 
homologous loci for its genes [19]. This evidently imposes a strict restriction limiting the 
availability of strings suitabh* for crossover. A second strategy proposed is to consider one of 
the parent string as the ’’master string” [22] and the other as ’’slave”. Before crossover the 
slave string genes are arranged in the same combination of the master string. The original 
combination of slave string is restored after the crossover operator is executed. Evaluation 
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of the fitness of new string evolved is then done by arranging the strings according to their 
locus. 


2.7 Differential Evolution 

Besides genetic algorithm a newer evolution technique has also been used as an optimiza- 
tion tool for obtaining the ground state structures of a-Si:H clusters. This method known as 
Differential evolution has been described by Kenneth Price and Rainer Storn [27]. In simple 
GA strings or the chromosomes (speaking biologically) are our starting material requiring 
a binary encoding. In DE the concept is different in the sense that real valued vectors are 
taken as the starting population and vector differentials instead of genetic differences is the 
driving force behind fitnc.ss improvement. An individual in DE thus comprises of an array of 
variables r('{)resented in real munbers. Starting with a random population the cost functions 
of the individuals are first evaluated. This step is followed by mutation, crossover and se- 
lection which are redc'fined in the presc'iit conte.xt. Any vector in the initial array is chosen 
as the parent. Two individuals are then randomly picked up from the given population 
with care taken to prevent selecting the same string twice. The difference between these 
two vectors are then calculated (see Fig. 2.3). Thus consider A' and Xi are the two selected 
vectors whose difference is giv(m as (A' - Xi). This difference vector is then multiplied with a 
parameter, known as DE mutation constant (F). The weighted difference vector thus created 
is added to a third vector Xj chosen randomly from the population forming another parent, 
which is given as 


A„,„ = .Y/+F.(V-Xi) (2.20) 

Crossover is then performed between vectors Xn^w the parent creating the trial vector. 
The mechanism of crossover involves substituting the variable in the parent with that of 
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Figure 2.3: Basic mechanism of differential evolution (Source: [27]) 


Xncw on the basis of some probability. A uniform random number is generated, which if 
greater than the crossover probability triggers variable swapping. This process is done by 
selecting the variables on a random basis obtained by performing a series of D-1 binomial 
experiments To ensure diversity in the population the last component to be changed is 
always taken from Xncw The trial vector thus evolved is then evaluated for its cost value. If 
the cost value is lesser than that of parent vector it is accepted into the next generation. This 
process is performed N times (X being the population size), each vector being selected once 
to play its role as a parent. Depending upon the problem to be optimized, the strategy for 
trial vector creation can be suitably changed. Further information in this regard is available 
in [27]. 


^For further elaborations see [27] 



Chapter 3 


Problem Formulation 


Alolecular structure optimization involves finding out the configuration for which total 
energy is minimum. Such a structure is energetically most stable and corresponds to the 
ground state structure. In order to achieve this through self-consistent calculations, the den- 
sity functional theory (DFT) by Kohn and Sham [2S] can be used, which however involves 
very large-scale computations. To avoid this computational overhead, efforts were made 
to calculate the total energy of clusters using semi-emi)irical models. Non-orthogonal tight 
binding theory is one of the commonly used model for simulating semiconductor systems like 
silicon. Here the total energy of the system is taken as : 

(R) = -b Uo (3.1) 

Where denotes the summation of repulsive pair potential terms Xij for the cluster and 
is given as: 


p^atom j_i 

E Exij 


1=2 j=l 


(3.2) 


The electronic energ}' is given by the expression 
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p^clec 

= J2 9k.ek (3.3) 

A:=l 

where Qk is the occupancy of the k’th eigenvalue e^-. Further details regarding the calculation 
of et and Xij 2 tre discussed in Sections 3.1 and 3.2. The correction factor Uq comprises of 
Uq’’ and Uq contributions for each of the atoms of silicon and hydrogen which is further 
described in Section 3.2. 

3.1 Calculation of the Eigenvalues 

According to the non-orthogonal tight binding theory the eigenvalue problem can be 
expressed as 


(3.4) 

Unlike in DFT the Hamiltonian in non-orthogonal tight binding theory (XOTBT) is not 
evaluated from the first principles. They arc obtained by using a model chosen on the basis 
of its ability to simulate the system. The eigenfunctions are then expressed in terms of a 
suitable basis set. Thus considering that |0') denotes the i'th function of the basis set the 
i’th eigenfunction can then be written as 

!'•■*> = ( 3 . 5 ) 

k 

It is implicitly assumed here that, the basis set ]i^*) chosen is sufficiently large and accurately 
represents the eigenfunctions of the system. The eigenvalue problem can thus be represented 
as 


^Ec“i?'") = £„Edr) 

m m 


(3.6) 
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Multiplying the above equation with {4>^\ we get 


m m 

(3.7) 

On defining the matrix elements Him and Sim as 


11 

(3.8) 

= (0'W") 

(3.9) 


the characteristic equation obtained is 

= 0 (3.10) 

m 

The above equation is solved for the eigenvalues using standard numerical techniques. 

3.2 Determination of Model Parameters 

Silicon atom consists of one s and three p orbitals while hydrogen has one s orbital. De- 
pending upon the type of orbitals and the nature of their overlap, four Hamiltonian matrix 
elements, B'sscr> f^pprr and four overlap integrals Ssp^, Sppa, Spp^: are defined. 

The matrix elements in non-orthogonal theory are taken according to [29]. If Vij denotes the 
Hamiltonian matrix elements in the orthogonal theory then Hij is given by, 

= ViAl + 1 - S|I (3,11) 

K 


If both atoms have s and p orbitals (e.g; matrix elements between two silicon atoms) then 
S 2 is given as 
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2^/^Ss 


^Sssa + 35 , 


PP7V 


(3,12) 


For one atom having s orbital and the other both s and p orbitals (e.g: matrix elements for 
a Si and H atoms) S2 is given by 


S2 = 



(3.13) 


When both atoms have s orbitals only (e.g: for two hydrogen atoms) , S2 is given by 
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The overlap integrals Sssa, Sspa, Spp^ and 5 pp,r arc given as 


(3.14) 




«(€a + fy) 


(3.15) 


where f aaV Hamiltonian matrix integrals in the orthogonal theory, k the non- 

orthogonality coefficient-efficient and e^, are the Hamiltonian matrix elements for the 
A and X' orbitals respectively. 

In this work, Hamiltonian matrix integrals of the orthogonal theory are assumed to have 
the form 


^ AA';.(c^o)ft^J^{-a(ry - do)} (3.16) 

where Vij is the distance between atoms and do the bond length parameter. The repulsive 
pair-potential term is expressed in the form : 


Xir-ij) = XQ€xp[-Aa{rij - do)] 


( 3 . 17 ) 
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The non-orthogonality coefficient-efficient k is given as 


~ ^0 + <^(^7 ~ do)'^ (3.18) 

where kq and a depends on the nature of atoms considered. The Hamiltonian matrix elements 
for the A orbital e\, depends only on the type of atom considered. Additional restrictions 
were imposed on the H-H interaction in the form of a cut-off function given as [18] 


/cut(^ 


1 + exp[ 


Xij - Tcut ' 
Scut 


(3.19) 


The off-diagonal terms in the Hamiltonian and H-H pair potentials were multiplied by this 
function [IS]. Further justifications for this can be found in [IS]. 

For silicon the Hamiltonian matrix elements ef', Cp* the bond-length parameter 

obtained from the work of Harrison [30] while the 
parameters S-re from the work of Menon and Subbaswarny [31]. 

The correction factor Uq^ for silicon has also been adapted from [31]. 

Parameters for Si-H interactions were obtained by fitting the experimental values of 
bond length, cohesive energy and vibration frequencies (A and E mode) for silane through 
least squares method [18]. This includes the parameters ef, 

and The H-H interaction parameters were obtained in a similar fashion 

using the experimental data on hydrogen molecule [18]. Parameters values used for each of 
the interactions are given in Tables 3.1 - 3.2. It is worth mentioning here that, the last term 
C/q in the total energy term is a correction factor, whose contribution in effect is taken as a 
constant throughout our calculation for the cluster. 
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Table 3.1: Interaction parameters for Si-Si, Si-H and H-H 


Parameters 

Si-Si 

Si-H 

H-H 

do (A) 

2.36 [30] 

1.48 [18] 

0.74 [IS] 

Xo (ev) 

0.05 [30] 

2.4365 [18] 

2.674 [18] 


1.70 [31] 

5.90 [18] 

3.30 [IS] 

a (A-') 

1.60 [31] 

0.7563 [18] 

1.200 [IS] 

^‘^sscr (^^0 

-2.37 [30] 

-4.02 [18] 

-13.250 [IS] 

^ spa ) 

2.52 [30] 

5.52 [18] 

- 

Vppa 

3.32 [30] 

- 

- 

^ppn 

-1.07 [30] 

I 

- 

'^cut (-'^) 

- 

- 

1.57 [18] 

^cut i*^) 

- 

- 

0.05 [18] 


Tabic 3.2; Parameters for Si and H atoms 


e. (ev) 

ep (ev) 

L'o (ev) 

-13.55 [30] 
-10.77 [18] 

-6.52 [30] 

1.000 [31] 
4.721 [IS] 


3.3 Application of GA 

Considering an atom to b(; a sphere with its mass concentrated at the center, any 
molecular structure is completely specified by the coordinate.s of the centers of its constituent 
atoms. In this work a cluster was considered to be situated within a cubic search space 
described by a Cartesian coordinate system. Thus in a cluster of type a total of 

m-f-n atoms can be explicitly described in terms of 3(m-l-n) coordinates. These 3(m+n) 
variables are mapped onto their corresponding binary forms, which on being assembled 
together gives rise to an individual ready for further genetic processing. Thus consider the 
following individual 


100101011101011 
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which represents the position of an atom in a cubic box of 125 cubic Avolume. The x, y, and 
z coordinates are represented by the 5 bit sub strings 10010, 10111 and 01011 respectively. 
Decoding the x, y, z sub strings into their real values using Eqn. 2.2, we get 

X = 2.90 A 
y = 3.71 A 
z = 1.77 A 

An individual is then checked for constraint violations prior to any further processing. The 
constraints imposed on a structure involve the minimum distances that must be maintained 
between Si-Si, Si-H and H-H atoms to prevent the atoms in the cluster from collapsing onto 
itself. Individual bond length parameters in Table 3.1 were generally used as the guideline. 
Attributing a low fitness value to an individual, which effectively rules out its transfer to the 
next generation, properly penalizes any violation of these constraints. Structures satisfying 
all these constraints arc evaluated for their colu'sive energy, the negative of which is assigned 
as the fitness. Ours being a minimization problem, an individual with the lowest cohesive 
energy receives the highest fitness value. 

Initially creating a random population generates a set of random structures whose indi- 
vidual fitness values are evaluated. Then the genetic operators discussed earlier are applied 
to this population giving birth to the individuals of the next generation. This process is 
continued till the population evolves to produce the fittest individual, which is then decoded 
to get the ground state structure. 



Chapter 4 


Computational Details 


4.1 Genetic Algorithm 

As discussed previously in Chapter 2, the search space in this problem consists of a cube. 
For all the calculations reported in this work, the side of cube was taken between 4 to 7 A. 
The string length for each variable was chosen such that the minimum resolvable distance 
remained around 0.01 A, or less. Thus for each variable, a 10-bit string was chosen at the 
least. Population size was kept in the range 90 - 100. Initially all the strings were niched for 
each of the variables using Goldberg’s triangular sharing function [19]. The objective was 
to prevent premature convergence and maintain an uniform selection pressure throughout 
the calculation as indicated earlier. This was followed by the application of reproduction 
operator. Tournament selection of size two coupled with elitist strategy was utilized for this 
purpose. The reproduction process was proceeded by the single point crossover operation. 
A crossover probability of 0.78 was generally used. Alternately, for bigger clusters (i.e for 
Siyli and SisH) two-point crossover was also applied with crossover probability in the range 
0.85 - 0.95. 

Crossover was followed by the mutation operation. Both jump and creep mutation 
operators w-ere exercised. Jump mutation probability was chosen in the range of 0.01 and 
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less. In all problems, the probability value for creep mutation was initially kept at 0.0011 
which was increased at higher generations. Depending upon the population variance at 
higher time steps, jump mutation probability was also adjusted to retain genetic diversity. 
For clusters in the series SimHn, total number of function evaluations required increases with 
the cluster size to be optimized. The overall strategy for optimization consisted of invoking 
runs from different initial conditions with the best result among all the runs being considered 
as the required global optimum. The stopping criteria for a run was kept at 1,00,000 function 
evaluations. 

For clusters SiH, SiH 2 , SiH^, SizH, SizHz, SizHz, SizHi, SizH and SuH, ground 
state results were obtained in less than 12,000 function evaluations. The optimum cohesive 
energy for clusters SiH.i and SizHz were deduced within 30,000 function evaluations. For 
clusters Si^H, SIqH and Si-H, a steep rise in the number of function evaluations were reg- 
isTt'red. Comparable results to simulated annealing [38] were achieved after approximately 
70,000 function evaluations. The SisH cluster posed the maximum difficulty with required 
function evaluations exceeding 1,00,000. This also indicates a non-linear increment in prob- 
lem stiffne.ss with increasing problem dimension. 

4.2 Differential Evolution 

Some of th(‘ clusters were also optimized using differential evolution algorithm. In this 
method tlu' s('arch space consich.'rcd was identical to that used in GA. Initially an assortment 
of random vectors was con.sidc'red. The population used was generally in the range 90 - 100. 
At the beginning of the search proce.ss, a high crossover rate (0.8 - 0.9) was applied in 
accordance -with the suggestion of Price and Storn [27]. This was followed by runs using a 
crossover probability in the range 0.5 - 0.6 with the F value varying between 0.6 - 0.7. As 
population converged, the crossover probability was reduced to a still lower magnitude of 
about 0.2 with a simultaneous decrease in the F value. The search process was continued 
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till the given stopping criteria was satisfied. 

All calculations in this work were done using a public domain FORTRAN 77 code [32] 
with several modifications. The suitability of the code had been previously confirmed in a 
separate study [33]. The entire computational job was performed in the local area network 
(LAN) at Indian Institute of Technolog}' Kanpur, predominantly using a Silicon Graphics 
server with central processing units : kllPS RlOOOO, Rev 2.6 under IRIX64 Released. 4 oper- 
ating system. The graphics output of the cluster geometries were generated using the XMol 
software on a SUN Solaris workstation. 



Chapter 5 


Results and Discussion 


Cohesive energies of the different ground state structures are given in Table 5.1 below. 
This table also provides a comparison with the previously published data in literature. 


Table 5.1; Ground state energy of various Si-H clusters 


Cluster 

Present Study 
(-ev) 

Simulated Annealing [18, 38] 
(-ev) 

SiH 

2.97 

2.97 

SiH2 

7.037 

7.04 

Silh 

8.90 

8.90 

SiH, 

13.409 

13.41 

SkH 

6.56 

6.57 

Si-Hh 

10.139 

10.14 

St2H, 

12.86 

12.87 

ShH, 

17.008 

16.93 

ShH-, 

17.907 

17.51 

Si-sH 

10.598 

10.624 

Si,H 

14.82 

14.83 

SkH 

19.385 

19.42 

Si,H 

23.895 

23.90 

Si-H 

28.139 

28.147 

SkH 

32.443 

32.445 


The bond lengths and bond angles of clusters are provided in Tables B.l - B.14. Struc- 
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tures obtained in this study are discussed in the following sections. In the figures given below 
the Si-H distances greater than 1.6 A but less than 2.0 A are shown in dark shades, while 
those still lesser are unshaded. Similarly for denoting Si-Si distance greater than 2.5 A and 
less than 3 A, dashed line has been used. The above choices are based upon the accepted 
covalent distances for these atom combinations. 

5.1 SiHi 


fnini 


1 



Figure 5.1: Ground state configuration of SiH 2 
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Silln exhibits a planar syininetrical structure. In this work the distance between the 

Si-II atoms is calculated as 1.51 A, and the ZH-Si-H is found to be 89.97°. These values are 
ill agreement with the experimental values of 1.5141 A, and 92° [34, 35] and the optimized 
results of simulated annealing [IS]. 

5.2 Silh 


1 



Figure 5.2: Ground state configuration of SiHz 


SiH-i possesse.s a pyramidal structure. Tiie distance between silicon and hydrogen resem- 
bles that for a covalent bonding with a bond length of ~ 1.5 A. The ZH-Si-H is approximately 
around 107° . 
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5.3 Si H 4 



Figure 5.3: GrouiKi state configuration of SiH.i 

SiH.i lias a fetrahetirai .struetun'. The optiniixed structure matches well with the struc- 
ture obtained I'.-qe-iim.'Uially [36j. This revalidates the parameters chosen for the purpose 
of our cohesive eiu'rgy calculation. Di.stance between the Si-H atoms is obtained as 1.48 A, 
and the calculated Zli-Si-H varies from 10G.4 to 111.6°. The cohesive energy values are well 
in acrord with the re.sult.s of [18] and [37]. 
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5.4 ShH 



Figure 5.4: (Jronnd state configuration of SiiH 

Si-^f! «‘.\iiib!{,s a iion-Hyuuiii'trieal structure. Si-Si distance is obtained as 2.34 A. The 
lone ii_vdr<j.!;en atriui i.s Ijonded strongly with one of the silicon atom, being situated at a 
distance of 1.54 A, which is similar to that obtained for covalent bonding. This structure is 
('(luivalcnt to the structure obtained by [ 18 ]. 
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5.5 SizI-h 



5.5: CJnnind stato ('onfiguratiou of Si 2 H 2 

In t!it* !\vn silirua atoms arc hoimci to each otiier with their bond lengths conipa- 

rablf to that found in Ihifh hydrogen atoms an^ bonded strongly to one of the silicon 

atoms with a bonding di.statietr of apjiroxi.’nafcly 1.54 A, and are oriented as mirror image 
to each other. I'he .stiuofnre oljtained is in agreement with the results of [18]. 
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5.6 Si-ol-h 



Figiu'i* ."i.fj: Ouniul state configuration of SUHs 

Till' >’/.>//} .structun* <it'tennin(Ml in this study shows that the bond length between 
its two Si atoms is greater than the accepted value for covalent bonding. This structure 
differs from Si, I!-, and Si-.H. l«)fh in the Si-Si bund length and also in the arrangement 
of iiydrogen atoms. Hydrogen atoin.s 4 and 5 arc situated symmetrically with respect to 
silicon atoms 1 and 2. Tiie tltird hydrogen atom is placed nearer to silicon atom 1 with the 
bond di.stance n semhliio. that of Si-H covalent bond. This structure closely resembles the 
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simulated annealing result of [IS]. 

5.7 Si2H.i 



6 


1 



I’ igui'c a.T; (^Iruund staler eonliguration of SltHi 

Stnulfire of >'?;//; eal-'ulated in tin- pre.stmt work Ls however, in variance with the sim- 
ulated annealiiii; results uf [IS]. ’The Si-Si distance and the arrangf'HKmt of hydrogen atoms 
closely re.sciu!»le to those in Siiih- While two hydrogen atoms are situated in symmetry 
wit is the silictm atoms, at a distance more than 1.7 A, the other two atoms are attached 
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to OIK' 
atoms 

5.8 


silicon atom each; their liistance being approximately 1.5 A. Both of these hydrogen 
are orii'utcd on the saiiK' side, unlike the structure reported by [18]. 

5 / ■>//:, 



Figure 5..S: CJrotmd state iMiifhoiratioa of Si2H5 


I'he cuiiesive t-nergy oluaiued for is -17.907 cv which is lower compared to -17.51 
cv obtained [FSh Tlie Si-Si bond length obtained is 2.995 A, compared to 2.29 A, reported 
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in the simulated aimoaling studies. Three hydrogen atoms are situated closely to one silicon 
atom, while the remaining two are attached to the other silicon. The Si-H bond length in 
each of tliese in.stanres is quite alike the accepted Si-H covalent bond length. 

5.9 



Figure 5.9; Ground state configuration of SizH 

In the SizH .structure .silicon atotn.s 1 and 2 are situated symmetrically \rith respect to 
tile tliird. The lone hydrogen atom is situated equidistantly from Si atoms 1 and 2. The 
Si-H bond distance being of the order of 1.8 A, which is greater than that existing in SiH 
cluster. Gupti* and Prasad [18] studied the characteristic frequency of the Si-H bond and 
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found its highost vibrational froqucnoy to be 1378 cm~^ 
in this duster is weak. 


This indicates that the Si-H bond 


5.10 SixH 


Ku’m 



1 


Figure 5.10; Ground state conHguration of Si^H 

The ground state structure of Si^H is quite similar to that of 5*4 duster [31]. The 
existence of hydrogen atom has only caused some distortion in the structure, retaining the 
essential geometrical feature of 5 * 4 . The hydrogen atom is situated at a distance of 1.54 A, 
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from the silicon atom 2 uhich is qt.itc alike that existing in Si-H covalent bond. Vibrational 
freqimncy study of th.’ Si-H bond by {3S] indicates the highest frequency as 2200 cm-\ which 
proves the pre.«ence of a strong Si-H bond for this .svstem. 


5.11 SiJI 



Figurt* 5.11: Ground .state configuration of Sifs,H 

ground .state structure as showti in the figure shows that the hydrogen atom is 
lying in the same plane as that of .silicon atoms 1,2 and 4. The hydrogen atom is situated 
chiscr to the silicon atoms 1 and 2 with Si-H distance being around 1.84 A. Vibrational 
spectrum studies (38j also indicate the absence of the characteristic frequency of strong Si- 



Figure 5.12: Ground state configuration of Si^H 

H biuid (21(10 arr^). The ground state structure is quite different from that found using 
simulated annealing [38]. It is evident that our results are a manifestation of a metastable 
state the eohesiv<? energy being 0.03 ev higher than that found by simulated annealing. 

5.12 SkH 

The structure of Si^H obtained in this work is similar to the results of simulated 
annealing [3S] and shows a remote correlation with the Sie configuration reported [39]. 
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Raghavadiari and Rohlfing [39] have suggested two alternate structures for Sie consisting of 
cither an edg(’-cap})ed trigonal bipyramid or a bicapped tetrahedron geometry. The structure 
of appears to bo a gro.ssl}- distorted bicapped tetrahedron, hydrogen causing the nodal 
silicon atom 3 to be displaced from its centered position in Sie. The optimized structure 
(ibtained through GA is 0.00-16 ev more energetic than that of simulated annealing. This en- 
ergy tlilh'ia'iict' is manifested in the variation of Si-H bond distances (± 5%). The maximum 
variation t'xists for the distance between atom numbers 2 and 7 with bond lengths differing 
by 0.10 A. However the Si-Si bond distances obtained are nearly same with the differences 
being lo.ss than ± 1% while the angular measures vary within ± 5 %. 

5.13 Si-H 

The optimized geometry for SirH predicted through GA differs from the results of 
simulated annealing [38] being 771 .Joules/molc less stable. Thus the result of GA indicates 
the existeuct' of a mettustable geometry for this cluster. The structure obtained however 
l)(iS.si'.s.M‘s ti gross similarity with the capped trigonal prism structure for Sir {C^v symmetry) 
considered by Rtighavachari and Rohlfing [39]. For the Sir structure thus proposed, the 
silicon atom 7 is arranged in a pyramidal fashion with atoms 1, 3, 5 and 6 constituting its 
l>a.sai plane (set! Fig. 5.13(a)). Tin? introduction of a hydrogen atom causes the atom 7 to 
shift towards one edge of the pyramid base. Comparing the ground state structure of SijH 
.^38'’ U) that of 5/7 (pentagonal bipyrarnid with Z? 2 /i symmetry) [31, 39], the hydrogen atom 
is fountl to position itself between silicon atom 7 and the mid plane comprising of 5 silicon 
atoms (see Fig. 5.i3(b)). It is to be noted that the optimized structure presented in this 
ivork difft'rs radictilly from the proposed ground state structure for ShH [38] in the absence 
of p<‘nta-bonded silicon atoms. 
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7 



(a) 


7 



(b) 


o.l,'?; (a) M('tastai)io coiifigurutioii of SijH obtained through GA (b) Ground State 
ctinligiiratiou of SirH ul)tuin(Hi through simulated annealing (Source: [38]) 

5.14 Si^II 

(b'uuiiil State l■oIl^i,',I;^afit)n for Si^H obtained through GA, conforms with the simulated 
aiinfaliiig results [38] with a cohesive energ}' difference of 0.0018 ev. The optimized structure 
b aftpr-i.vimately similar to that of Sin- .Menon and Subbaswamy [31] has proposed a distorted 
bicas.pe<! octalu dron structure as the lowest energ}' state for Sk. The addition of a hydrogen 
auuu uiudifii‘H this geometry further. Thus the H atom seems to push silicon atoms 6 and 
8 (.see Fig. 5.M) furtlier apart, compared to their positions in Sk. The hydrogen atom 
tiiereby, lii's nearest to atoms 6 and 8 at a distance of ~ 1.75 A which is much greater than 
the Jifcepted covalent distance for Si-H. The next nearest neighbors i.e., atoms 1, 2, 3 and 4 
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5.14; Ground state configuration of S'UH 
are piacft! at a distance of ~ 3.0 A. 

5.15 Results of DE 

Diffejentia! evolution (see section 2.7) method has been used &zi alternate evolutionary 
optinii/ation teciiniqiie in this work. Initially the smaller clusters of SimHn series have been 
optinii/ed a.s a test ca.S(>. The cohesive energies of clusters thus obtained, are listed in Table 
5.2. Frelimitiary results arc comparable to those of GA. Further work with this approach is 
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reciuirt'd to establish the relative advantages or disadvantages of DE over GA. 


Table 5.2: Results of Differential Evolution 


Cluster 

Results of 
DE(ev) 

Results of 
GA(ev) 

SkH 

6.56 

6.57 

SkHi 

10.09 

10.139 

ShHz 

12.85 

12.87 

SiH2 

7.04 

7.04 

SiH^ 

8.90 

8.90 






Chapter 6 
Conclusions 


Ground state structures of SinHm clusters for (n= 1,2 and m= 1,5) and SiiH clusters 
for (1=3,8) have been optimized using a tight binding non-orthogonal scheme coupled with 
genetic algorithm. Most of the results obtained in this study compare well with those 
obtained from available experimental work and the optimized results of simulated annealing 
[18. 38). The structures obtained for Si^Hs, Si^H and SirH&ve the only ones, which 
.«:igtufi<‘untly differ from the work of Gupte and Prasad [18, 38]. The cohesive energy values 
dotermitied for and are lower than those obtained through simulated annealing, 
remlering the structures reported in this work thermodynamically more favorable. Using 
a .similar line of retusoning, structures of SkH and SijH are expected to be metastable 
coinpartfd to those predicted by Gupte and Prasad [38]. 

The SitH scries of structures displaj's some special features. It is found that the hy- 
drogen atom in each of them resides outside the framework of silicon atoms being nearly 
cquidi-stantly positioned from its closest neighbors. However any interaction between them 
i.s expected to be weak since the mutual distances are much greater than that of Si-H cova- 
lent bond. This was conclusively proved in the vibrational frequency spectrum studies [18]. 
These structures are also found to be quite sensitive to small energy changes. Thus a energy 
difference of 0.008 ev is enough to change the geometry drastically from the ground state of 
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Further stiulics on these clusters and their derivatives are required to improve our un- 
derstanding about such emerging materials. Thus using the ground state results, molec- 
ular dynaniiCvS simulation can be conducted to get an inkling about the various physical 
and cheniieal proi)erties of these clusters. It is expected to facilitate the understanding of 
hydrogen-induced effects in a-Si:H as mentioned previously. This work also establishes the 
effect iveiK’ss of GA in obtaining the ground state structures with much reduced computa- 
tional efforts. This can be extended further by parallelizing the code for genetic algorithm. 
A successful venture is expected to open newer vistas of research, especially in the fields of 
materials simulation and nanostructures a subject of immense importance for the future. 
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Appendix A 

Coordinates of Different Si-H Clusters 


Table A.l: SiH 2 


Alum"' 

X-Coordinale 

Y-Coordinate 

Z-Coordinate 

Sii 

2.9402 

3.0300 

0.2400 


2.2201 

1.7099 

0.0510 

H, \ 

4.3197 

2.3700 

0.1500 


Table A.2: ShH 



X-Coordiuatc? 

Y-Coordiiiate 

Z-Coordinate 

■Vn 

1 

1.4471 

0.1653 

2.2582 

1.7649 

1.1057 

0.1188 

H, 

1.8940 ‘ 

2.4590 

0.8470 


Table A.3: SiH^ 


Atom 

X- Co ordinate 

Y-Coordinate 

Z-Coordinate 

Sii 

1.6280 

2.4889 

1.5420 

H-z 

0.3051 

1.9960 

1.0330 

Hz 

1.8550 

3.8429 

0.9470 

H, 

2.7008 

1.5891 

0.9940 


Table A.4: Sm 


Atom 

X-Coordinate 

Y-Coordinate 

Z-Coordinate 

6’/. I 

1.4212 

1.3712 

1.3710 

Hz 

0.4821 

2.3718 

0.7910 

Hz 

2.2479 

2.0380 

2.4090 

//, 

2.3717 

0.8520 

0.3460 

Ih 

0.6550 

0.2590 

1.9760 


Table A.5: SkHi 


Alum 

; X-Coordinate 

Y-Coordinate 

Z-Coordinate 

Sil 

2.5018 

1.4980 

3.1750 

Siz 

2.5490 

1.4980 

0.8560 

Hz 

1.4981 

2.6630 

2.9870 

\ //■, 

1.3491 

0.4730 

2.9990 


Table A.6: Si^Hz 


dum 

X-C3)or<Unate 

Y-Coordinate 

Z-Coordinate 

'Sii 

1.4580 

0.3341 

1.7790 

S'U 

1.7909 

2.7299 

0.2840 

hI 

2.7051 

0.0000 

0 .9760 

H.i 

1.9760 

2.0010 

1.9640 

* *1 

0.7910 

1.1611 

0.4199 






Table A.7: S^Hi 


Atom 

X-Coordinate 

Y-Coordinate 

Z-Coordinate 

Sh 

3.4S30 

1.4381 

1.8690 

Sh 

4.7360 

3.4441 

3.4541 

Hz 


4.3838 

2.2700 

H, 


2.7790 


Hs 


1.8201 


He, 

3.7569 

2.3190 



Table A.8: SkH-^ 



X-Coordiiiate 

Y-Coordinate 

Z-Coordinate 

Sn 

! 1.0759 

1.3008 

1.7907 

Si 2 

3.4540 

2.9353 

2.5930 

Hz 

2.6711 

1.2432 

1.8201 

Hz 

3.7380 

4.0698 

1.6140 

Hr, 

2.9160 

3.9238 

3.6400 

He, 

0.S219 

1.9082 

3.1411 

[Jk 

0.7240 

2..3390 

0.7342 


Table A.9: Si^H 


’ X-Courcliaate 

Y-Coordinate 

Z-Coordinate 

1.9589 

Sis 

4.5510 

2.2459 

2.7508 

1.8360 

0.8394 

1.1204 

2.8041 

3.0969 1 

3.0030 

3.4129 









Table A. 10: S^H 


Atom 

X-Coordinate 

Y-Coordinatc 

Z-Coordinate 

Six 

3.2940 

4.6274 

2.7843 

Sk 

3.9805 

2.4705 

3.4510 

Si, \ 

1.07S3 

4.7059 

2.0392 

Si, i 

1.8627 

2.5296 

2.4118 


4.7647 

2.1373 

2.1765 


Table A. 11: Si^H 


.Atom 

X-Coordinate 

Y-Coordinate 

Z-Coordinate 

Sti 

3.2314 

3.5137 

1.9610 

Sh 

3.1370 

0.8785 

1.5216 

Si, 

1.3491 

2.5096 

1.0664 

! 5... 

1.2863 

2.1802 

3.4197 

I 5 .’-, 

3.6707 

1.9451 

3.6235 

1 //,’ 

1.7.569 

0.4866 

2.6510 


Table A.12: Si^H 


Vtum 

X-C.’oordinate 

Y'Coordinate 

Z-Coordinate 

"sir 

1.0165 

1.3637 

2.8690 

Si'j 

4.4770 

2.6540 

1.5836 

Si, 

1.S817 

3.5875 

2.8104 

Si, 

2.5367 

1.3294 

1.0215 

Si, 

3.31S7 

1.4467 

3.4409 

c, 

OUt 1 

4.0909 

3.6950 

3.7047 


2.4976 

3.3089 

1.1143 



Table A.13: SijH 


Atom 

X-Coordinate 

Y-Coordinate 

Z-Co ordinate 

Si, 

4.9506 

3.7401 

40978 

Sl2 

2.9487 

3.7431 

5.4373 

Siz 

1.3131 

3.1979 

3.4148 

Sii 


5.2555 

4.6048 

Sh 


5.3112 

3.3297 

Si, 


1.6795 

4.1973 

Si- 

4.0948 

2.5032 

2.2599 

ll, 

2.6937 

3.8339 

2.1690 


itoni f X-Ctjonii 


Si 2 
Si -2 
Si, 
St-, 

S 1 1 ', 
Sir 

C!'* ' 
if |h, 

//., 


1 .(M3' 
5.4781 
5.2901 
2.938 
4.01 r 


,3548 


4.1761 




Appendix B 


Structural Details of Si-H Clusters 


Table B.l; SiHi 


Bund 

Bond la'iigth 

Angle 

Bond Angle 

Six ~ //■> 

1.5144 

H 2 - Six - Hz 

92.62 

Six - H, 

1.5319 

Hz -Hz- Six 

44.00 

Hi - //:, 

2.2027 

Hz - Hz - Six 

43.38 


Table B.2; SiHz 



Bund Length 

Angle 

Bond Angle 

s' ) . 

* ' 1 

Hi 

1.5007 

Hz - Six - Hz 

107.22 

: Six 

- Hz 

1.4962 

Hz - Six - Ha 

107.98 

Six 

- Hi 

1.5036 

Hz - Sii - Hi 

106.76 

Hi 

- Hz 

2.4125 

Hz- Hi- Sii 

36.33 

Hi 

- Hi 

2.4303 

Hi -Hi- Sii 

36.04 

[Jk 

- Hi \ 

2.4078 

Hi -Hz- Sii 

36.45 
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Tabic B.8; SioHs 


Bond Length 


2.9951 

1.5965 

3.8451 

3.6994 

1.5023 

1.5225 

2.0184 

1.5251 

1.537C 

2.S78C 

3.3561 


Angle 
Si2 - Sii - 
Si2 - Sii - 
Si2 - Six - Hs 
Si2 - Six - He 
Hr - Six - Si2 
He - Six - H, 
He - Six - He 
He - Six - He 
He - Six - Hr 
Hi - Sk - He 
Ih - SU - He 


Bond Angle 
38^87 
21.50 
23.68 
70.90 
89.84 
48.32 
61.28 
99.61 
105.56 
116.74 
91.38 


Table B.9: ShH 


St I Sie 
Stx Ih 
St 2 •' Six 
St, ~ Hi 
St. - Ih 


Humi L<'ngih 
2.775': 
2.3G3‘J 
1.7761 
2.360f 
1.809.' 
3.00U 


Angle 


Sh - Six - Sie 
Sh - Six - Hi 
Sie -Six -Hi 
Six - Sk - Hi 
Sie - Sh — Hi 
S i ^ **” S te "— Hi 


Bond Angle 
53.98 
39.72 

91.85 

38.85 
91.09 
36.25 










Table B.13: SijH 


Bond 

Bond Length 

Angle 

Bond Angle 

Sii — Si2 

2.40SS 

Sk — Si, — Sis 

45.05 

Sh - Sh 

3.7406 

Sk — Sil — Sii 

33.65 

Sil — Sii 

3.9784 

Sk — Sil — Sk 

66.95 

Si, - Sis 

2.3984 

Sk — Sil — Sk 

53.79 

Six - Sk 

2.8442 

Sk — Sil — Sk 

97.55 

Si, ~ Siy 

2.3749 

Sk - Sil - Hs 

74.31 

Si, - Ih 

2.9702 

Si^, — bii *5 ’^4 

35.70 

Sl 2 — Sis 

2.65761 

Sk — Sil — Sk 

51.02 

Si 2 — Sii 

2.3821 

Sk — Si, — Sk 

39.81 

S >2 - Sis 

2.6516 

Sk — Si, — Sir 

55.44 

.,- 1 , At • 

bi'2 - oui 

2.4078 

Sk — Sil — Hs 

31.47 

Sii — Sir 

3.5981 

Sii - Si, - Sk 

33.34 

Sii - //« 

3.2794 

Sii — Si, — Sk 

68.96 

Si,, — Sii 

2.3770 

Sii — Sil — Sir 

88.10 

Si., - Sk 

2.9082 

Sk - Si, - Hs 

51.31 

Si.i -- 5/o 

2.3956 

S i^ ~~ S ii — Siq 

90.84 

Si , Sir 

3.0909 

Si:i — Sil — Sir 

81.18 

Si., -Ik 

1.9653 

Si;, — Sil — Hs 

41.52 

Sii - Sir> 

2.3742 

Sk - Sil - Sir 

53.24 

Si 1 Si a 

3.9742 

Sk — Sil — Hs 

61.48 

1 Si, '■ Sir 

4.5655 

Sir - Sil - Hs 

40.55 

Si 1 - Ik 

3.1426 

Sil — Sk — Sk 

95.05 

Sf, - Sk 

3.7473 

Si, - Sk - Sk 

112.28 

Si:. Sir 

3.1059 

Si, — Sk — Sk 

56.34 

Siu - Ik 

1.0770 

Sil - Sk — Sk 

72.39 

Si,, - Sir 

2.3763 

Sil - Sk - Sir 

40.87 

l (j 1 1 ^ 

2.9740 

Sil - Sk - Hs 

60.69 

' Si: - Ik 

1.9344 

Sk — Sk — Sk 

55.96 



Table B.14; SkH 


Hoiui 

Roiul Length 

Angle 

Bond Angle 

Sii - SI 2 

4.0275 

Si 2 — Sil — Siz 

34.74 

Sii- Sfi 

3.903S 

Sio — Sil — 8 ( 4 . 

40.96 

Sii - su 

2.3704 

S 12 S^i S 25 

34.79 

Si\ - Si-, 

3.7855 

S 12 S^i Szg 

35.47 

Sii - Si, 

3.8000 

Si: - Sil — Sij 

41.91 

Sii — Sij 

2.3985 

Si: — Sil — Sis 

55.00 

Sii - Sin 

2.3836 

Si:- Sil - Hq 

48.82 

1 Sii - //.j 

3.2579 

Siz ““ Sil — Sii 

75.70 

1 Si. - Si:i 

2.3706 

Siz — Sil — Sis 

59.15 

i Si. — Sii 

2.7241 

Siz — Sil — Sis 

59.04 

1 Si-: - Sir, 

2.3472 

Siz - Sil - Sij 

35.58 

' Si-: - Si, 

2.3942 

Siz ““ Sil — Siz 

35.01 

1 — 517 

2.7560 

Siz — Sil — f/g 

53.19 

1 5;, ~ 5;s 

3.3001 

Si.i -Sil- Sis 

36.41 

1 Si. - //., 

3.0912 

Sii - Sil - Sis 

37.35 

, Si,i -■ Sii 

4.0356 

Si.i -Sil- Sir 

70.86 

i Si.i - Sir, 

3.7966 

Si.j ““ Sil ”** Si'Z 

88.00 

1 Si.i - Si, 

3.7969 

Si-i — Sil — Hd 

64.34 

j Si.i - Si: 

2.4004 

Sis -Sil- Sis 

60.50 

' Sti ■■■ Si, 

2.3829 

Sir, — Sil — Sij 

37.76 

Si-i ~ III 

3.2578 

Sis — Sil — Sis 

88.45 

\ Sii ■ Stu 

2.3461 

Sis - Sil - Ho 

81.78 

1 Sii ■ Si, 

2.3955 

Sk - Sil - Sij 

77.32 

Sl\ Si-; 

2.7646 

Sis — Sil — Sis 

54.45 

1 Si, 

3.3023 

Sis - Sil - Ho 

27.01 

i Sil - ' II, 

3.0894 

Sh - Sil — Sis 

70.53 

i 

: Si:, Stf, 

3.8212 

Sh - Sil - Ho 

84.05 

1 Si-, Si; 

2,3930 

Sis - Sil - Ho 

32.08 

: Si’, ■' Si, 

4.4185 

Sil — Si: — Siz 

69.78 

Si-, II, 

4.6281 

Sh - Si: - Sil 

34.78 

'■ Si,-, Si; 

4.0240 

Sil — Si: Sis 

66.96 

' Sh. ■■■ Si.* 

3.0965 

Sil - Si: - Sis 

67.07 


St,-. - 11, 

1.7306 

Sh - Si: - Sh 

35.54 



2.7609 

Sil — Si: — Sis 

36.28 


Sn ih 

3.8400 

Sil — Si: — Ho 

52.49 


Si* ~ Ih 

1.7708 

Siz - Sh - SU 

104.56 
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